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Abstract 

In neutrino oscillation with non-standard interactions (NSI) the system is enriched with CP 
violation caused by phases due to NSI in addition to the standard lepton Kobayashi-Maskawa 
phase 5. In this paper we show that it is possible to disentangle the two CP violating effects by 
measurement of muon neutrino appearance by a near-far two detector setting in neutrino factory 
experiments. Prior to the quantitative analysis we investigate in detail the various features of the 
neutrino oscillations with NSI, but under the assumption that only one of the NSI elements, 
or Eer; is present. They include synergy between the near and the far detectors, the characteristic 
differences between the and Eer systems, and in particular, the parameter degeneracy. Finally, 
we use a concrete setting with the muon energy of 50 GeV and magnetized iron detectors at two 
baselines, one at L = 3000 km and the other at L = 7000 km, each having a fiducial mass of 50 
kton to study the discovery potential of NSI and its CP violation effects. We demonstrate, by 
assuming 4 x lO^i useful muon decays for both polarities, that one can identify non-standard CP 
violation down to |ee/i| — ^ few x 10 and \£eT 

I ~ 10-2 at 3(7 CL for 6I13 down to sin^ 26*13 = lO""^ 
in most of the region of 5. The impact of the existence of NSI on the measurement of 6 and the 
mass hierarchy is also worked out. 
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I. INTRODUCTION 



Discovery and establishment of neutrino mass and the lepton flavor mixing [1] by the 
atmospheric, the solar, and the reactor experiments triggers a tantalizing question of 
whether neutrinos have interactions outside the Standard Model of particle physics. In fact, 
the possibility has been extensively discussed starting from the early works in jl, 0, [s], 0] 
and in conjunction with expectation of new physics at TeV energy scale 0, [!]■ 

To our current understanding, one of the most signiflcant features of fundamental matter 
is the quark-lepton parallelism. Then, the crucial question is to what extent it prevails and 
in what respect it fails For example, quark and lepton flavor mixing is described by 
the CKM and the MNS [l| matrices, respectively, similar in structure but with very 

different values of the mixing angles [l2|. It is well known that in the quark sector the 
Kobayashi-Maskawa (KM) theory lll| for CP violation works perfectly to the accuracy of 
the experiments achieved to date |13|. We still do not know if a natural extension of the 
KM theory to the lepton sector provides a correct description of nature (l4j |. 

In the presence of non-standard interactions (NSI) of neutrinos, however, the task of 
exploring lepton CP violation is inevitably confused or enriched by the coexistence of possible 
another source of CP violation induced by NSI [l5|. In this paper we discuss the problem 
of how to resolve the confusion caused by the two sources of CP violation in the context 
of measurement of the lepton mixing parameters in neutrino factory experiments [igI]. Of 
course, discovery of NSI will open an entirely new window to physics beyond the Standard 
Model of particle physics. Therefore, in this paper, we aim at achieving the two goals: 
Illuminating how to resolve the confusion between the standard interaction (SI) and the 
NSI parameters, and exploring the discovery potential of NSI. In this paper, we mean by 
SI parameters and 6 assuming that the other standard lepton mixing parameters are 
determined to reasonable accuracies. 

The problem of confusion between the SI and the NSI parameters has been raised in the 



form of ^13— NSI confusion in |17l. |18[|. In a previous paper [19j we have shown that it can be 
resolved by a two-detector setting, one at baseline L = 3000 km and the other at L = 7000 
km, in neutrino factory experiments, enjoying an intense neutrino beam from a muon storage 
ring. In a natural continuation of the work, we discuss the problem of two-phase confusion 
in this paper by extending our treatment to allow the NSI elements 8^/3 (for deflnition, see 
Sec. ini) to have phases. For a related work on the same subject see [13]. The similar question 
of distinguishing two kind of CP violation in the context of "unitarity violation" approach 



2l| has also been investigated [2^, [23 



One of the most important questions on NSI is how large their effects are. Given the 
extreme success of the Standard Model at low energies the natural framework to address the 
question is via higher dimensional operators [13]. Assuming the new physics scale of M^p ~1 
TeV, it may be ~ {Mw/M^pf ~ lO'^ and ~ (Mw/Mnp)"^ ~ 10"^ for dimension six 
and eight operators, respectively. Possible forms of theses operators are further narrowed 
down by the severe constraints on four charged lepton processes [25|, the SU{2) partners of 
the neutrino interactions on leptons, the problem raised in [8|. In view of the fact that the 
operators which are free from the four lepton counterpart are very limited [1^, it may be 
the right attitude to anticipate and cover a wide range of magnitudes 10""^ < leapl < 10~^ 
(or less if possible) as our target for hunting the NSI effects. For this reason, we try to cover 
the region down to Sajj ~ 10~^ in our analysis in this paper. The constraints on magnitude 
of NSI on neutrinos has been investigated [23, 28|, and the authors of 29] made major 
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advances in understanding the loop constraints. 

As is well known, the effects of NSI exist not only in pro pag ation in matter but also 
in production and detection processes of neutrinos 0, 11, 18, 30 1. In this paper, however, 
we only deal with the effect of NSI in neutrino propagation in matter. The reason for this 
limitation is mainly technical; The features of neutrino flavor transformation are sufficiently 
complicated so that it deserves separate treatment. For the same reason, we restrict ourselves 
to systems with a single NSI element, either e^fj, or Ser, at one time. Thanks to these 
limitations, we can have a transparent view of the parameter degeneracy [3l|, 32, 33 



m 

systems with NSI [3^ in its limited aspect that can be seen more clearly under the simplified 
setting. In this paper, we will try to give a complete understanding of the characteristic 
features of the system, not just showing the sensitivity contours. In particular, we will clarify 
the reasons why the sensitivities are so different between the systems with and 

Taking j£p,^ I ~ 10~^ — 10~^ as the target region, it is unlikely that reactor 
perbeam [36| experiments are powerful enough to have significant discovery potential for 
NSI. Thus, we are left with either the neutrino factory [l6| or the beta beam [33| as the 
possible technologies within our current knowledge. In fact, there exist a number of articles 
which are devoted to discuss the capabilities of discovering or constraining NSI in a neutrino 



1351 or su- 



factory experiment, for instance, see Ref. [17|, Il8|, Il9|, l30|, l38 



Because of the clean event 
reconstruction capability, we focus on the appearance channels z/e — z/^ and Ue — ^ the 
so called golden channels for the neutrino factory [3^. We assume, following [19], a high 
energy muon beam of E = 50 GeV and two magnetized iron detectors at L = 3000 and 7000 
km. Notice that in the i/g-induced appearance channels, this or z/g u^-, the only relevant 
NSI elements are either eg^i or as shown in a perturbative framework in jl^], which is 
nothing but the natural extension of the one developed in [39|. 

Of course, it is important to place constraints on NSI from all possible available means. 
In fact, the bounds (to be) placed on NSI have been discussed in the context of accelera- 
tor neutrinos [13] (excluding the neutrino factory), atmospheric neutrinos |4l|, reactor or 
spallation beam neutrinos [i^], solar neutrinos [43j, as well as of astrophysical neutrinos IIZ 



II. FEATURES OF NEUTRINO OSCILLATION WITH NSI 

Neutrino oscillation in the presence of NSI is a highly nontrivial problem. We try to give 
an introductory discussions on it in the light of some recent progresses. 



A. Neutrino propagation with NSI; General framework 

We consider NSI involving neutrinos of the type 

= -2V25gG^(z7„7^P^z.^) (77^'P/), (1) 

where is the Fermi constant, / stands for the index running over fermion species in the 
earth, f = e,u, d, and P stands for a projection operator which is either = 1(1 — 75) or 
= 1(1 -|- 75), and a, P = e, fi and r. 

To summarize its effects on neutrino propagation it is customary to introduce the effective 
e parameters, which are defined as SafB = X]/p^^«^5 where n/ (ng) denotes the /-type 
fermion (electron) number density along the neutrino trajectory in the earth. Using such e 
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parameters the evolution equation of neutrinos in the flavor basis is given by 
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where U is the MNS matrix [l] (for which we use the standard notation [l^]), o, = 2\/2G pn^E 
j^, is the neutrino energy and Amf - = — m| with nii (2 = 1 — 3) the neutrino mass. 

Here are some cautionary remarks: We note that since our discussion in this paper ignores 
effects of NSI in production and detection processes of neutrinos, the results here must be 
interpreted with care. Yet, the system is complicated enough so that the analysis of NSI 
effects in neutrino propagation itself seems to be a meaningful step. We also note that 
since we concentrate on NSI effects in neutrino propagation in matter only the vector type 
interactions can be probed. Moreover, the types of NSI that can be studied may be limited 
to a sub class of generic NSI [45| . 



B. Neutrino oscillation probability with NSI 



We focus on z/^ 



z/^ and Ve 



appearance channels to analyze effects of NSI, anticipat- 



ing intense z/g and Ve beam from either a neutrino factory or beta beam. All the quantitative 
analyses will be done assuming muon storage ring and magnetized iron detectors with muon 
charge identification capability. For simplicity, we do not consider disappearance channels 
in this paper as we believe that the inclusion of them would not modify substantially our 
results. For qualitative understanding of various features of the analysis we use the oscilla- 
tion probability derived in e perturbation theory js^]- We assume that 6'i3 is small as it i s 
natural for the neutrino factory setting, though the current bound is still rather mild 
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In this framework, e = 



Am: 



Am: 



sin 6 
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keal (tt = A*, t) are regarded as small expansion 
parameters of the same order, while is regarded as of order unity. 

It is shown that to second order in e the relevant NSI elements are only eg/x and e^r in 
the Ve channel (a = /x, r). Then, the oscillation probability in the 

can be written in the form of an absolute square as 



z/^ channel 
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where = cos 6ij , Sij = sin 6ij , 



and A31 = 



As we will learn in the next subsection, the system with Sf,^ and e^T is too complicated 
to do a full analysis, having six unknown parameters including three CP violating phases. 
In this paper, therefore, we work on systems with only a single type of NSI, e^^ or Eer, as a 
step to understand the features of neutrino oscillations with NSI. 
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The expression of the oscillation probability in ([3]) illuminates a notable difference between 
the systems with Ee/j. and Eer- Namely, the sensitivity to Ee/j. is generally better than that to 
Est because the two kinematic factors in the square bracket multiplied by the latter (former) 
tend to cancel (add up). 



C. Bi-probability plot in neutrino oscillation with NSI 

Unfortunately, the reduced system with a single type of e is still quite complicated. 
Therefore, we try to understand the characteristic features of neutrino oscillations with E^a 
(a = /i, r) in this subsection. What makes this system so complicated is the presence of two 
CP violating phases, the lepton KM phase [ll| in the MNS matrix [l[ , and the phase 0ea of 
the NSI element Eea = l^eaW'^"" (« = An appropriate tool to illuminate the properties 

of neutrino oscillation with emphasis on the role played by phases is the bi-probability plot 
— P(z/e — > Vy) (for short P — P) space proposed in H . 



The appearance oscillation probability in the neutrino channel can be generically written 

as, 

P = P{ue ^ Uf,) = A + {71 + Bcos6 + C sin 6) cos + ( J - i3 sin 5 + C cos 6) sin 0, (4) 

where the expressions of the coefficients in (jl]) can be easily derived from ([3]), and are given 
in Appendix |X1 The corresponding expression for anti-neutrino channel can be obtained by 
doing the transformations a —>■ —a, 6 —>■ —6, and (p ^ In principle, there are two ways 
to draw the bi-probability plot, varying (p while holding 6 fixed, or varying 6 while holding </> 
fixed, the "traditional" way. (See flAip in Appendix [Al) The expression in (jl]) anticipate the 
former way. We use both ways for convenience to help to illuminate the physical properties 
of the system. 

Let us start to proceed step by step to understand the system with NSI. In Figs. [1] and [21 
we present the bi-probability plots for the systems with Eefj, or Eer, respectively. We use the 
former way, varying while holding 6 fixed, to draw these figures. The neutrino energy 
is arbitrarily chosen as E = 20 GeV, but we have checked that the main features of the 
bi-probability plots are similar for other values of the energy in the region 10 GeV < E < 
40 GeV, which are relevant in our setting. The values of NSI parameters in both systems 
are chosen so that the appearance probabilities take values comparable with each other. 

There are some distinctive features in the bi-probability plots which are notable in Figs. [T] 
and [2] including difference between the systems with Ee^ and Eer- 

• The slope of the major axis is mostly positive and negative in the Eer and e^^ systems, 
respectively, apart from the region around 6 = ir. Furthermore, ellipses rotate when 
S is varied. It is the very origin of the dynamical behavior of the system with two 
CP violating phases which spans a wide range in the bi-probability space. The reason 
for this behavior is explained in Appendix [Bl and this property will be important 
to understand the difference between the Eer and Ee^ systems in sensitivities to the 
parameters. In particular, when the two phases are varied, the ellipses span over 
almost the whole triangular region in P — P space in the Eefj. system. See the left panel 
of Fig. E] in SeclYl 

• A glance over Figs. [1] and [2] indicates a global feature, the size of the ellipses strongly 
(mildly) depends on \e\ (sin^ 26*13). Also, the relative distance between ellipses is 
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Bi-Probability plots for e^^ , L=3000km, E=20 GeV for sinlQ^^ = 10 \nd lo'^ 



o 



I I I I I I I I I I I I I I I I I I I 
'o - 



le 1 = 0.0005 

eji.' 



^1 I I I I I I I I I I I I I I I I I 



t 

a: 

I >, 



1 1 1 1 1 1 1 1 1 1 1 1 
- le 1 = 0.0005 


1 1 1 1 1 1 1 


:% H - 

1 1 1 1 1 1 1 1 1 1 1 1 


5 = 

5 = 71/2 J 

5= 71 

5 = 371/2 - 
1 1 1 1 1 1 1 



I I I I I I I I I I I I I I I I I I I 
: le 1=0.001 O. 



^ I I I I I I I I I I I I I I I I I 



■•©I 3 
2 



1 1 1 1 1 1 1 1 1 1 1 1 1 
: le 1=0.001 


1 1 1 1 1 1 






o. 




— cf^. 








1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 



I I I I I I I I I I I 



le 1 = 0.005 




I I I I I I I I I I I I I I I 

le 1 = 0.005 




I I I I I I I 



12340 1234 

P(v -v^) [xlO"'] 

FIG. 1: Bi-probability plots for the systems with £e^ drawn by continuously varying (j^e^^ fo^ four 
different values 5 = 0,7r/2,7r and 37r/2. The panels in the upper and lower rows correspond to 
sin^ 2^13 = 10~^ and 10~^, respectively, whereas the ones in the left, middle and right columns 
correspond to |ee^i| = 5 x 10~^, 10~^ and 5 x 10~^, respectively. The normal hierarchy was assumed 
and the values of the other mixing parameters used are sin^ = 0.31, Am^i = 8 x 10^^ eV^, 



sin O22, = 0.5 and | Am|-|^| = 2.5 x 10 eV . We use these four values for the rest of the paper. 
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Bi-Probability plots for e^^ , L=3000km, E=20 GeV for sin^2ej3 = 10 ^ and lO"'* 
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FIG. 2: Similar plots as shown in Fig. [1] but for the systems with Se,-. The panels on the left, 
middle and right columns correspond to jeerl = 5 x 10"'^, 10~^ and 2 x 10^^, respectively. 
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greater for smaller \e\ and larger sin^2^i3. In regions where overlapping between 
ellipses is significant we must expect stronger two-phase confusion and parameter 
degeneracy. Whereas in regions where ellipses are sparse, we expect less degeneracy 
and better sensitivities for the determination of parameters. 

We will revisit these features of the bi-probability plot in our overview of the sensitivities 
in Sec. llVt thereby meriting understanding some characteristic features. 



D. Near-far two-detector setting in a neutrino factory 



In this subsection, we briefly review the two-detector setting in a neutrino factory for 



determination of the NSI and SI parameters proposed in [1^. We denote, hereafter, the 
detectors at baselines L ~ 3000 km and L ~ 7000 km the near and the far detectors, 
respectively. The basic idea behind it is that the far detector, which is located close to the 
so called "magic baseline",^ plays a role complementary to the near detector;^ The synergy 
between the two detectors greatly strengthens the potential of parameter determination and 
help solving the confusion between the NSI and the SI parameters. This is very similar to and 
concordant with the idea of using the far detector as degeneracy solver in the measurement 
of the standard oscillation parameters without NSI 3l|, |47| ■ The setting is now considered 
to be a "standard one" in an international effort for designing a neutrino factory [55| . 

In our setting of turning on a single NSI element, Eefi or Est, the oscillation probability 
has special features at the magic baseline. It takes the form in the system with Sefi of 



^■'23 ■'13 



a — ^mfi 



sm 



+ 



cos(5 + 0e^) + S23a\e 



e/i I 



sm 



AE 



(5) 



P(z/e — ^ z/^; Est) cau bc obtained by the transformation 523^6^* C23SeT at the magic base- 
line, in accord with the prescription given in Appendix El Notice that the two phases come 
together in the form S + (pefi in ([S]); The physical phase must be unique at the magic base- 
line because the solar Am^i effect is absent and the system is effectively two generations.^ 



The naming is due to 47|. The distance has been known as the matter refraction length ^3|], as can be 
seen, for example in 48[ . The authors of [49[ pointed out that 6 (CP phase) dependence goes away at the 



magic baseline. One of our motivations for placing the second detector at the distance in 19| was that it 



is the best place, roughly speaking, to detect the effect of matter density to which the NSI elements are 
proportional [sol. 

^ Though proposed in a variety of contexts, the basic idea for the two-detector setting here is quite different 



from the one in reactor ^13 experiment 35|, in which the cancellation of systematic errors between the 
detectors is the key element. Ours here is to seek a complementary role, and hence it is closer to the one 
for measurement of CP violation 51]. The idea has also been applied to the Kamioka-Korea identical 
two-detector complex to determine the mass hierarchy as well as discovering CP violation (a possible 
option for upgrading the T2K experiment 53) in which these two aspects are unified 



53 



54|. 



The reasoning was spelled out in [19], and it has subsequently obtained an analytic proof as a general 
theorem of phase reduction in generic systems with all NSI elements j34| . 
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Clearly, the bi-probability plot shrinks into a line with only a cosine dependence on the 
CP phase, similar to the case at the effective oscillation maximum [s^], and the length of 
the shrunk ellipse is proportional to \eea\- This behavior is clearly seen in Fig. [3l and its 
consequence in some of the allowed contours in Sec. HV] for example, in Fig. [7| and Fig. [HI 

Here, let us add a comment; The baseline L = 7000 km chosen is not exactly the magic 
baseline which is more like L = 7200 km. We have chosen this distance anticipating that 
it is not always possible to have a site which coincides to the exact magic baseline. The 
difference between the two baselines is illustrated in Fig. [31 

Bi-Probability plots for L = 7200 and 7000 km for sin^2e = 10"^ 
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FIG. 3: Similar bi-probability plots as shown in Figs. [J and [2] (where (j) is varied continuously) 
but at the magic baseline -Lmagic = 7200 km (upper panels) and at the actual baseline used in this 
paper L = 7000 km (lower panels), for sin^ 2^13 = 10~^. The left (right) panels are for the case 
where |ee/j| = 0.001 (|eer| = 0.001). Comparison between the upper and the lower panels indicates 
the difference of the behavior at Lmaeic = 7200 km and L = 7000 km. 





III. ANALYSIS METHOD 



We describe the method which will be used in our analyses in the following sections. In 
this paper, we make the same assumptions as we made in our previous paper 1^ for the 



neutrino factory set up. We assume an intense muon storage ring which can deliver 10^^ 
useful decaying muons per year with the muon energy taken to be 50 GeV, two magnetized 
iron detectors with a fiducial mass of 50 kton each, one at baseline L = 3000 km and the 
other at L = 7000 km. In this work, we assume that each detector can receive neutrino flux 
corresponding to 10^^ useful decaying muons per year. We also assume data will be taken 4 
years in neutrino mode and another 4 years in anti-neutrino one. We consider only the golden 
channels, z/g —>■ and z/g in this paper. For simplicity, we use the constant matter 
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density approximation, and take the earth matter density along the neutrino trajectory as 
p = 3.6 g/cm^ and p = 4.5 g/cm^ for baselines L = 3000 km and L = 7000 km, respectively. 
The electron fraction Ye is taken to be 0.5. We believe that using more realistic earth matter 
density profile will not change much our results. 

We always use the normal mass hierarchy as an input for the results shown in this paper, 
unless stated otherwise. We, however, vary over the two different mass hierarchies during the 
fit of the data. We will consider the setting in which only one element of NSI, either e = e^ij, 
oi e = Eeri where e stands for \e\ and its complex phase 0, is turned on. For simplicity, we 
fix the values of oscillation parameters relevant for solar and atmospheric neutrinos (apart 
from the mass hierarchy which is assumed to be unknown) to their current best fit values as 
sin2^i2 = 0.31, Amii = 8 x 10"^ eV^ sin^ ^23 = 0.5 and jAm^^l = 2.5 x 10"^ eV^. The last 
approximation may affect the estimated sensitivities quantitatively, but not in a significant 
way. 

Our function, which is similar to the one used in [l9|] apart from the part which include 
systematic uncertainties and background, is given by, 

>f ^ V ^^^''^ ~ ^.S°(^i3. 5. sign(AmlJ)] ' 

where N°j^i^ is the number of observed (simulated) events computed by using the given input 
parameters and N^^j^ is the theoretically expected number of events to be varied in the fit by 
freely varying the standard mixing and NSI parameters. The summations over the indices 
i,j and k imply summing over 3 energy bins, 2 baselines (3000 km or 7000 km), and 2 modes 
(neutrino or anti-neutrino), respectively. The intervals of 3 energy bins considered are 4-8 
GeV, 8-20 GeV, and 20-50 GeV for neutrinos and 4-15 GeV, 15-25 GeV, and 25-50 GeV 
for anti-neutrinos. We refer the readers to Ref. 1^ for details on how and A^/j^fc are 

computed. 

Using the function defined in Eq. (Q, we define the allowed regions in terms of 2 
parameter space (shown in Sec. irv]) by the commonly used condition, Ax^ = ~ xLin 
= 2.3, 6.18 and 11.83 for 1, 2 and 3 a confidence level (CL), respectively, for 2 degrees of 
freedom (DOF). For the sensitivity plots which will be shown in the next section, we used 
the condition Ax^ = ~ Xmin = 1)4 and 9 for 1, 2 and 3 a CL, respectively, for 1 DOF. 

To illuminate the physical properties of the systems with NSI and to indicate discovery 
potential of the two-detector setting in neutrino factory we think it is revealing to utilize 
the following two choices for the experiment related assumptions in the function in ([6]): 

Choice R: We take, following ||55i], cXsys = 2.5 % for the signal and ctbg = 20 % for the back- 
ground systematic uncertainties. We assume an energy independent detection efficiency of 
70 % for both neutrino/antineutrino modes. This is a rough approximation to the efficiency 
curve obtained in [H^, which is presented in the left panel of Fig. 20 of this reference. The 
background fractions (efficiencies) are also taken from the same reference, presented in the 
right two panels of Fig. 20 of (56| . 



Choice O: As a second choice we assume no systematic uncertainties, no background, 100% 
detection efficiency. Namely, N^j^i^ = and CTgys = ctbg = 0. Only the statistical errors for 
signal events are implemented. Vanishing systematic errors is nothing but an approximation 
for a wide range of situations in which they are much smaller than statistical uncertainties. 

The symbol "R" is meant to be "realistic" , while the "O" optimal. 
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We use the results of the analyses with two different recipes of the experimental uncer- 
tainties to make our presentation most informative. In Sees. [IV] and |Vl we rely on the 
results of analysis based on Choice O of the parameters, that is, no errors except for the 
statistical one. In these sections we try to illuminate the structure of neutrino oscillation 
with NSI in the golden channel. Hence, we prefer to deal with the allowed regions which are 
free from obstruction by the experimental errors. 

Whereas in Sees. |Vl] and IVIII where the sensitivity contours of the NSI and the SI pa- 
rameters are discussed, we present the contours for both choices. Choice R and Choice O of 
the parameters. At the minimum, the comparison between these two choices reveals the 
effects of the systematic uncertainties and the efficiencies on the sensitivities, and therefore 
is useful. If we are to prepare for search for effects induced by dimension eight operators 
beyond dimension six ones, as argued relevant in Sec. [T] (see also [Hsj), one may be interested 
in optimal sensitivities that can be achieved by an ultimate detector.^ While the concrete 
design of such a detector is not known, its first approximation may be given by taking the 
limit of vanishing errors under the current detector setting. 

IV. OVERVIEW OF THE SENSITIVITIES AND SYNERGY BETWEEN THE 
TWO DETECTORS 

In this section we give an overview of the sensitivities to the NSI elements, |£:eQ,| and 
(j)ea (a = /^^t), and the SI parameters, 5 and sin^2^i3, to be achieved by the detectors at 
L = 3000 km and L = 7000 km separately and in combination. The principal purpose of 
the discussions in this section is to understand the global features of the sensitivities of the 
complex systems. In fact, various viewpoints have to be involved to really understand what 
the complementarity between the near and far detectors means: 

• How prominent is the synergy between the near and the far detectors for determination 
of SI and NSI parameters? How it differs between the systems with Se/j. and ^er? What 
about the dependence on values of the parameters, in particular on the size of the NSI 
elements? 

• How can the two-phase confusion be resolved? Does the answer depend on which NSI 
element is turned on? 

• What is the nature of the parameter degeneracy in the system with NSI, and whether 
it can be resolved by the two-detector setting? If so how can it be realized? 

The first two points will be discussed in this section and we will devote the whole Sec. IVlfor 
the problem of degeneracy. 

A. Data set 

Toward the above stated goal, we systematically generated the contours of allowed regions 
determined by the data taken in the near (3000 km) and the far (7000 km) detectors and the 



^ We give a caution here that the current estimation of experimental systematic uncertainties and back- 
grounds may not be sufficiently mature to address the question of sensitivities to such tiny values of 
NSI. 
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one combined based on the analysis method described in Sec. IIIII We call the set of allowed 
contour figures the "data set" . To generate the data set we take the input parameters in the 
following way: sin^ 2^13 = 10"^ and IQ-^, 6 = 0, 7i/2, n and 37r/2, = 5 x 10-^ lO'^ and 
5 X 10"^, 0eAi = 7r/4, 37r/4, 5tt/4 and 77r/4. Altogether, there are 2 x 4 x 3 x 4 = 96 sheets.^ 
Then, the above parameters are varied to fit the data. During the fit, unless otherwise 
stated, the parameter space of both the normal and the inverted neutrino mass hierarchies 
are swept. The input is always taken as the normal hierarchy except for the results shown 
in Sec. IVIIBI The allowed regions at 2a and 3a CL are defined by the region corresponding 
to "~ Xmin < 6.18 and 11.8, respectively. We try to illuminate the dominant features of 
the complex system of neutrino oscillation with NSI by analyzing the data set. 

In the system with Ser, because of poorer sensitivities compared to the system with Sg/i, 
we use an order of magnitude larger values of the NSI parameters to generate the data set: 
l^g^l = 5 X 10"^, 10~^ and 2 x 10~^, (per = vr/4, 37r/4, 57r/4 and 77r/4, with the same values 
of 6^13 and 6. The ranges of the input values of and Eer are chosen such that they lead to 
the similar values of the probabilities as can be seen by comparing the bi-probability plots 
of both systems in Sec. Ill CI They roughly correspond to the sensitivity regions obtained in 
"l9|. Again there are 96 sheets as in the Eefi case. Due to the lack of space, we only show 



a tiny fraction of such allowed contours in the following subsections but below we will give 
quahtative summary of our data set described above. 



B. Classification scheme of the data set 

By reviewing 2 x 96 sheets we recognize that there exist 4 typical types: 

• Type A: Measurement at the near detector (3000 km) cannot produce closed contours 
but the combination with the far detector (7000 km) makes the contours closed (at 
least for one of the 4 variables), which determines the parameters within errors. This 
is the case where the synergy is most prominent. 

• Type B: Measurement at the near detector does produce closed contours, but multiple 
solutions are allowed, this is the phenomenon of parameter degeneracy. In most cases 
in this category the combination with the far detector solves the degeneracy. 

• Type C: Measurement at the near detector itself produces closed contours without 
degeneracy. The role of the far detector is to make the allowed region smaller. 

• Type D: Both the measurement at the near detector and the combination with the far 
detector fail to produce closed contours. 

The closing or not of contours, of course, depends upon the confidence level. We use 2 and 
3 a CL in the tables and figures in this section. 

In Table HI and HI] we present number of sheets which fall into the categories of Type A— D 
of the data set prepared for the system with NSI element Ee/j, and Eer, respectively. In the 



^ A sheet, in principle, contains six two-dimensional plots because there are four parameters to be fit. 
Practically, however, three of them appear to be enough; We have used the two-dimensional plots with 
axes |£e^| - 4>ef_i, 6 - 4>ei_i, and sin^ 26*13 - i^- 
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type and magnitude of NSI 


sin^ 26*13 


Type A 


Type B 


TypeC 


Type D 


\ee^,\ = 5 X 10"'' 


10-4 


5(2) 


5(8) 


1(5) 


5(1) 


lee^l = 5 X 10"^ 


10-3 


3(2) 


5(1) 


3(8) 


5(5) 


ke^l = 10-3 


10-4 


0(0) 


15(8) 


1(8) 


0(0) 


ke^l = 10-3 


10-3 


2(4) 


4(3) 


6(8) 


4(1) 


lee^l = 5 X 10-3 


10-4 


4(0) 


12(14) 


0(2) 


0(0) 


|ee^| = 5 X 10-3 


10-3 


0(0) 


11(11) 


5(5) 


0(0) 



TABLE I: Presented is the number of sheets which fah into the categories of Type A— D out of 
96 sheets of the data set prepared for the system with NSI element Eg^. The definition of Type 
A— D is described in the text. The numbers in and out of parentheses are the ones at 2 and 3(T 
CL, respectively. The results shown here (or in this and next sections) correspond to the choice O 
(100% detection efficiencies, no background and no systematic uncertainties) defined in Sec. IIIII 



type and magnitude of NSI 


sin2 26*13 


Type A 


Type B 


TypeC 


Type D 


= 5 X 10-3 


10-4 


16(16) 


0(0) 


0(0) 


0(0) 


\£er\ = 5 X 10-3 


10-3 


16(16) 


0(0) 


0(0) 


0(0) 


\eer\ = 10-2 


10-4 


16(15) 


0(1) 


0(0) 


0(0) 


\£er\ = 10-2 


10-3 


14(11) 


1(4) 


1(1) 


0(0) 


|eer| = 2 X 10-2 


10-4 


16(15) 


0(0) 


0(1) 


0(0) 


\eer \ = 2 X 10-2 


10-3 


8(4) 


6(9) 


2(3) 


0(0) 



TABLE II: Same as in Table H] but for the system with NSI element e, 



system with Ser the classification scheme works except that there is no Type D as we will 
see below. Here, we need some details of the criteria for the classification. Existence of a 
closed contour^ does not necessarily mean that the unique solution is contained in it; It can 
occur that multiple solutions are involved in a single region because they are too close to 
be resolved due to limited statistics, for example. Our classification scheme A— D above is, 
therefore, purely "phenomenological" in nature. 

It should also be noticed that classification of figures into A and D has an ambiguity; It 
occurs in some cases the allowed contour is closed for \e\ but not for sin2 2^i3, for example. 
We define Type D such that the allowed contours of any one of the four variables fail to be 
closed after combining the yields from the two detectors. Consistently, Type A is defined 
such that at least one of the four variables has closed contours. It makes classification scheme 
a complete one; The number of pages add up to 16 when summed over Type A— D for a 
given input values of and sin2 26*13. 

A glance over the tables illuminates the following features: 

• In the system with e^^- most of the sheets are classified into Type A except for the 



^ In doing computation, for technical reasons, the regions of parameter scan were limited to IQ-*^ < \se^\ < 
6 X IQ-'^ and 5 x IQ-"' < |eeT-| < 3 x IQ-^ in respective systems. Therefore, closed contours imply they 
are closed in the above regions. 
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case of relatively large ^13 and leerh Synergy between the two detectors are far more 
prominent compared to the system with Se^. 

• In the system with Se^ synergy effect exists but its effectiveness depends on 61^ and 
\Se^\ in an intricate way. There is a general tendency that closed allowed region is 
more common in cases with large \Se^\- For a given value of \Sep\ synergy is stronger 
in case of smaller 6*13, however with the notable exception of the case \se^\ = 10~^. 

One may ask why the synergy is so powerful in the Eer system. The answer involves two 
reasonings; Measurement at the near detector is not as powerful as in the system because 
the two kinematic factors in the coefficient of e^t in ([2D, which have similar magnitudes, tend 
to cancel, resulting the lower sensitivity to Eer- Then, powerfulness of the synergy is due to 
the higher statistics of the far (7000 km) detector due to larger {Eerl, which is enhanced by 
the former reasoning. 



C. Features of sensitivities in system with eg/^ 

In this and the next subsections we give further description of the features of sensitivities 
to the NSI elements, \Eea\ and (pea (« = A^, t), and 6 and sin^ 26*13 by the near-far two-detector 
setting. In what follows we present only one figure for each category A— D. The category 
B (the case with parameter degeneracy) will be discussed in the Sec. |Vl We use 3cr CL 
criterion for the classification of types. 

Let us first consider the system with Egfj.- As an example of Type A, we present in Fig. H] 
the regions allowed by measurement at the near detector [L = 3000 km, top panels), the far 
detector (L = 7000 km, middle panels), and the two detector combined (bottom panels). 
It is for the input values of the NSI element \Ee^\ = 5 x 10~^ and (pe^ = 77r/4 whereas the 
SI parameters are taken as 5 = vr and sin^ 26*13 = 10^"^ with the normal mass hierarchy. 
Left and right panels in this figures show allowed regions plotted in \Eei_i\ — 4>en and 6 — 0e/i 
space, respectively. We note that for this case, no allowed regions exist for the inverted mass 
hierarchy regime. 

In this figure, it is indicated that measurement at the far detector, though it itself does 
not appear to be so powerful, helps closing the allowed regions, indicating prominent synergy 
between the two detectors. This feature is consistent with the one observed in systems with 
NSI but with frozen phase degree of freedom discussed in the previous paper [l9| . However, 
Type A is relatively minor in the data set for large values of \Ee^\ as seen in Table [H indicating 
less prominent synergy between the two detectors in the Eefj. system compared to the EerS. 
In Fig. m accuracy of determination of 0e/i is limited because of small \Ee^\, though 6 — (pe^ 
confusion is resolved. 

An example of Type C is presented in Fig. [5l which was generated using the same values 
of \Ee^\ and 6*13 as in Fig. H] but different phase values, 5 = 37r/2 and (pefi = 57r/4. For 
relatively large values of \Eef_i\, the case of closed contours as the result of near detector 
measurement is most popular, as can be seen in Table [11 In fact, as long as NSI parameters 
are concerned, there is no example, within the cases examined, of unclosed contours (only 
upper bound) at 3000 km detector for the input value |ee^| = 5 x 10~^. In the example in 
Fig. [3, we still see modest synergy between the detectors. Interestingly, the sensitivity to S 
improves due to the far detector measurement at the magic baseline. The S — 0e^t confusion 
is clearly resolved though the accuracy of determination of 0e/i is niodest. 
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FIG. 4: An example of Type A. Allowed regions in the 



■'e/j. 



'6/1 



plane (left panels) and 



plane (right panels) corresponding to 2 and 3 a CL obtained for the system with eg/i- Upper 
(middle) panels correspond to the case where only a single detector at 3000 km (7000 km) is 
taken into account, whereas the lower panels correspond to the case where the results from the 
two detectors are combined. The input parameters are taken as: sin^ 2^13 = 10~^, 6 = tt, \£eii\ = 
5 X 10"'^ and cpe^ = Yir/A (indicated by the green asterisk), and the input mass hierarchy is normal. 
No allowed regions exist in the inverted mass hierarchy regime. All the results shown here (or in 
this and next sections) correspond to the choice O (100% detection efficiencies, no background and 
no systematic uncertainties) defined in Sec. Illli 



In Fig. [6] presented is an example of Type D generated by using the same parameters as 
used in Fig. [5] except for (pen = 77r/4, showing that the synergy is least prominent among the 
categories A— D.^ In this case the unclosed contours at 3000 km fail to be closed at 2a CL 
even after combination with 7000 km detector. Again the far detector improves sensitivity 
to 6, but not at all to (pen- 

Finally, it may be worth to note that the characteristic feature at the magic baseline that 
only the combination 4>efi + ^ is constrained [l^ is barely visible in most of the data set. The 
possible reasons are: The baseline is not exactly equal to the magic baseline. The statistics 
at the far detector is of course much less and sensitivities to the parameters are close to the 
sensitivity limit, lacking clear indication of 6 — (p correlation. However, appearance of 6 — (p 
correlated oblique strip becomes frequent in case of relatively large \een\ = 5 x 10~^. 



^ We note that for sin^ 2^13 = 10 ^ (not shown), some cases classified as Type D at 3a CL become Type 
A or B even at 2a CL. 
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FIG. 5: An example of Type C. Same as for Fig. [5] but for the input parameters sin^ 2^13 = 10~^, 
6 = 37r/2, |ee^| = 5 x 10~^ and i^e/i = 57r/4, and the normal mass hierarchy. No allowed regions 
exist in the inverted mass hierarchy regime. 



D. Feature of sensitivities in system with Ser 

Synergy between the two detectors is much more 



er than the system. To demonstrate this point we present 

n/^nc im I cr I A, 

-1-3 



We now turn to the system with e 
prominent in the system with e 

examples of Type A in Fig. [7] the allowed regions in Iserl — (per and 6 — (per space for the 
input parameters, l^erl = 5 x 10~^, (per = 77r/4, sin^ 2^^i3 = 10^^, 6 = 3n/2 and the normal 
mass hierarchy. Within the currently used parameters most of the data set with Ser fall 
into Type A except for the largest values of |£er| and 613 in the set, l^erl = 2 x 10~^ and 
sin^ 2^13 = 10~^, as indicated in Table HTl We note that in the system with Ser, for the case 
where only the near detector is used, the clone solution in the inverted mass hierarchy (the 
input is the normal one) exists very often, as shown in the second row of Fig. [3 This will 
significantly lower the sensitivity to the mass hierarchy determination compared to the case 
where no NSI of the e^r type is present, which will be demonstrated in Sec. IVIIBI Type C 
is rare in the data set in the system with for the currently used values of the parameters, 
as indicated in Table HTl An example is shown in Fig. [8] which was generated with the input 
parameters Iserl = 2 x 10~^, (per = 77r/4, sin^ 2^13 = 10"'^, 6 = 11/2 and the normal mass 
hierarchy. 

Synergy is very remarkable in the system with Ser apart from the case of the largest values 
of l^erl and 6'i3 mentioned above, as clearly seen in Table HTl Yet, effect of the far detector 
data is visible even in relatively large \eer\ case, making allowed regions tiny as a result of 
combination of both detectors as shown in Fig. [HI This feature is generic for larger values 
of \eer\- The allowed contours in the (per — S plane in Figs. [7|and[H]due to the far detector 
measurement alone clearly displays characteristic feature that only the combination (per + S 
is constrained. 

It is remarkable to find that there is no Type D in the data set, representing an extremely 
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powerful synergy between the two detectors in the system with Ser- It is also very significant 
that 6 — (per confusion is resolved in most of the figures for the system with Ser- 

Though we have achieved reasonably good understanding of the features of the sensitiv- 
ities, there is a curious feature in Tables [T] and [TTl Let us compare Fig. [T] and Fig. [2] in the 
right-bottom panel with sin^ 2^13 = 10~^, = 5 x 10~^, and £er = 2 x 10^^. The values of 
the NSI parameters are chosen so that the size of the probabilities are comparable. Then, 
one would expect that the (relative) accuracy to which we can determine (or constrain) 
NSI parameters is similar between the Eefi and the Est systems; The ellipses with different 
6 overlaps and the degree of overlapping is very similar in both systems. Nevertheless, the 
potential for parameter determination is in fact very different between Sen and Est systems 
as one observes by comparing between Table [T] and Table [Til Type B is dominant in the Ee/i 
system, whereas Type A dominates leaving no Type B in the E^r system. It clearly indicates 
that the sensitivity is far better in the ^e^i system compared to the Eer system. On the other 
hand, the parameter degeneracy is much severer at the near detector in systems with Eefj, 
compared to the ones with shown in the tables. 

A unified understanding of the puzzling features becomes possible once one draws the 
bi-probability plots by varying 6 for four (twenty) different values of (^13) as done in 
Fig. m The degeneracy is severer in the Eefi system because of the more dynamic behavior 
of the bi-probability ellipses as shown in the left panel of Fig. [91 Because the ellipses can 
locate themselves essentially everywhere in the bi-probability space there are chances that 
fake solutions can be produced at points far apart from the true solution. What about the 
difference in the sensitivities? We observe in the right panel of Fig. [HI that the ellipses in the 
Eer system remain in a much more compact region when S and 6*13 are varied. Because of 
the finite resolution of the experimental data it appears that the dense concentration of the 
ellipses with different parameters leads to merging of many degenerate solutions, resulting 
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FIG. 6: An example of Type D. Same as for Fig. [Dbut for the input parameters sin^ 2^13 = 10 
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FIG. 7: An example of Type A. Allowed regions obtained for the system with Ser- The input 



parameters are taken as: sin 26 



13 



10 



-3 



5 X 10-'^, 6 = 3tt/2, and 



7tt/4 and normal 



mass hierarchy. Unlike the plots shown in Figs 01161 for the case where only the single detector at 
3000 km is taken into account, allowed regions exist not only for the normal mass hierarchy regime 
but also for the inverted one, as shown in the panels in the second row. 



in the lack of the sensitivity in the Eer system. 

By comparing Figs.lHElfor (which show no strong synergy by combining two detectors) 
and Figs. [71 and [H for Ser (which show strong synergy), one might think that strong synergy 
for Eer is simply because of the larger values we considered for this NSI element than that for 
Eefi- We note, however, that if we consider larger value for Egfj., the allowed NSI parameters 
obtained from the near detector alone are already restricted to small regions and we do not 
see any strong synergy when combined (see Fig. [T2l) with the far detector unlike the case of 



V. PARAMETER DEGENERACY WITH NSI 



In this section we try to understand the property of the degeneracy in systems with NSI, 
encountered in the previous section, i.e., in the category T ype B. It was indicated that 
the structure of the conventional parameter degeneracy jsil. Is^. \3^ . (for reviews see e.g.. 



|49|, |59j), the intrinsic, the sign-Amg^, and the 6*23 octant degeneracies prevails but in a 
way involving NSI parameters, in particular, the phases of e parameters [SJ]. In fact, the 
degeneracy in systems with NSI turns out to be a highly complicated problem. In a general 
two Eea (« = Z^, t) system there are total six parameters. Though some progress has been 
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FIG. 8: An example of Type C. Allowed regions obtained for the system with Ser- The input 
parameters are taken as: sin^ 2^13 = 10~^, leerl = 2 x 10~^, 6 = 7t/2, and (j)eT = 77r/4 and normal 
mass hierarchy. No allowed regions exist in the inverted mass hierarchy regime. 



made in 3J], its complete treatment still eludes us.^ In our simplified system of a single 
type of e, however, we can achieve a good understanding of the phenomenon. After giving 
some examples to show characteristic features in this section, we give in Appendix O some 
semi-analytic treatment to obtain the degenerate solutions for a given true solution. 



A. 0-degeneracy 

We start from the simplest case which may be called as the ^-degeneracy, in which 
the degenerate solutions are characterized only by distinct values of 0, and values of the 
remaining variables are nearly equal. In this subsection we restrict ourselves to the system 
with Sen- In Fig. [To] some examples of such degeneracy are presented. 

Let us understand why and how such degeneracy arises. We start from the expressions of 
appearance oscillation probabilities in (jlj). Because the terms proportional to C and I (and 
corresponding coefficients for anti-neutrino, C and X) are an order of magnitude smaller than 
the other terms for typical values of the parameters for relevant for neutrino factory, they 
can be neglected. Then, it is interesting to examine a simplified model for the oscillation 
probabilities 

P-A= ±v^^B^T^^TZB^TOsIcos(0 + a) 

P-A= ±Vb^ + n'^ + 2BTZ cos 5 cos(0 + a), (7) 



A completely new type of degeneracy, the solar-atmospheric variables exchange degeneracy, is found 
there. But, this degeneracy does not appear to enter into the discussion in this paper, and hence it is 
outside of the scope of this paper. 
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FIG. 9: Bi-probability plots drawn by continuously varying 5, for four different values of (j) for 
L = 3000 km. For each value of (j), the 20 different values of sin 26113 between 10"^ and lO'^ 
(logarithmically varied) are considered. 
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FIG. 10: Some examples of (/)-degeneracy. The normal mass hierarchy is assumed. The parameters 
are taken as sin^ 2^13 = 10~^, \£eij\ = 5 x 10~'^ and the baseline is L = 3000 km. 



where the upper (lower) line corresponds to neutrino (anti-neutrino) channel and the sign 
must be chosen correctly, a and a are given by 

BsmS _ Bsm5 , , 

tana = — — — , tana = -= ;: — (8) 

Scos^ + 7^' -Bcos5 + 7e ^ ' 

The features of the simplified model ([7]) can be best illustrated by the bi-probability 
plot in P — P space. See the top-right panel in the upper figure in Fig. [T] in particular, 
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the case of S = 0, for one of the clearest example. Namely, the ellipse shrinks into line 
due to the presence of only cosine dependence of the phase. Then, it is evident that the 
degenerate solution exists; If we have a solution at = 0i then we have another solution at 
02 + « = 27r — (01 + a), namely, 02 = 27r — 0i — 2a. 

The simplest situation arises at 5 = where a = a = 0. The second solution is given by 
02 = 27r — 01 which nicely explains features of the left panel of Fig. [TDl In the right panel 
of Fig. [ini the input value of 6 is and a can be obtained by using ([H]) as a = —0.88, 
where a typical value of the energy E = 30 GeV is used. Then, corresponding to solution 
01 = = 2.36 there is a second solution 02 = 27r — 0i — 2a; = 5.7, which explains the right 
panel of Fig. [10] in a reasonable accuracy. 

Though this type of degeneracy can be solved in most cases with the far detector mea- 
surement, it can fail in certain cases. Since dependence at the magic baseline ([5]) takes a 
particular form cos (5 + 0) the 0-degeneracy cannot be solved for 5^0. An example of this 
phenomenon is given in the left panel of Fig. [TOl 



B. Intrinsic and Am|i-sign flipped degeneracies 

Now, we discuss more generic form of the degeneracy, the intrinsic and the Arn|i-sign 



flipped degeneracies which are similar to those in the standard neutrino oscillation [3l|, l32 
The 023 octant degeneracy [33| does not arise because we take the maximal ^^23- 



In Fig. [TT], we show, for the system with Se^, the case (which is categorized as Type B) 
where the measurement at the near detector at 3000 km leaves four regions at 2 a, two of 
which correspond to the normal hierarchy and the other two the inverted ones. It is the 
manifestation of the sign-Amgi degeneracy in systems with NSI where the phase of NSI 
parameter is indeed actively involved. Unlike in most cases this is one of the special 
cases in which the far detector measurement fails to resolve the intrinsic degeneracy. In this 
particular case it occurs in the following way; Because the relation cos((52 + 02) = cos((5i + 0i) 
approximately holds the far detector cannot solve the degeneracy. Also, the energy spectra 
are very similar for the two solutions with the same sign of Am'^^. 

In Fig. [12] an example of Type B for the system with with degenerate solutions at 
the near detector is presented. As in the system with Type B takes place when l^erl 
is relatively large. Here, again the sign-Am|i degeneracy exists at the near detector (top 
panels), but it is completely resolved when combined with the far detector measurement. 
Most likely, the intrinsic degeneracy does not exist in this case, or is already resolved by the 
near detector measurement. 



C. Does spectrum information solve the degeneracies? 

One may ask the question, "to what extent are the degeneracies we saw in the previous 
subsection robust?" It is well known that the intrinsic degeneracy in the conventional 
neutrino oscillation with SI is often fragile to the spectrum analysis. To give an insight to 
this point we present in Fig. [12] the oscillation probabilities as a function of neutrino energy 
which correspond to the same Amgi-sign intrinsic degeneracy (left panel) and the flipped 
Amgi-sign degeneracy (right panel). From these figures it is obvious that the degeneracies 
we observe in the systems with or Ser are quite robust against the spectrum analysis 
as long as we consider only a single detector at L = 3000 km. One can go through the 



20 



10" 



10" 



10'; 

10" 



1 1 1 1 1 1 I 1 1 I 1 1 I 1 1 I 1 1 I 
* 

2a 

3a 



normal 

I I I I I I I I I I I I I I I I I I I 



inverted 

I I I I I I I I I I I I I I I I I I I 



2.0 
1.5 
1.0 
0.5 
0.0 
1.5 
1.0 



.1 I I I I I I I I I I I I I I I I I I. 
_ normal _ 



: i I I I I I I I I I I I I I I i : 



inverted - 



10" 



10" 




I I I 1 1 1 1| 1 — I I I 1 1 1 1| 1 — n: 

normal - 



I I lllll| — I I I lllll| — h-H 
inverted ■ 



I I I I L_i; 



sin^2ej3= 10 ^ 5 = 71/2 
le | = 5xlO"\(|) =71/4 

I I I I I I I I I I I I I I I I I I I 



— * — 
1 1 / 1 1 1 1 1 1 1 1 1 1 1 1 


.1 1 M M| 1 1 l_ 


- 1 1 1 1 1- 



0.0 0.5 1.0 1.5 

d) Ik 



10 



10^ 



sin^2e, 



FIG. 11: An example of Type B. Allowed regions in the (j)efj, — l^e^il plane (left column), 
plane (middle column) and sin^ 2^13 — 6 plane (right column) corresponding to 2 and 3 a CL 
obtained for the system with Se^- Panels in the upper 2 rows (3rd row) correspond to the case 
where only a single detector at 3000 km (7000 km) is taken into account, whereas the ones in the 
4th row correspond to the case where results from the two detectors are combined. The input 
parameters are taken as: sin^ 2^13 = 10~^, 6 = 7r/2, \ee^\ = 5 x 10~^ and (/>e^ = '/r/4 (indicated 
by the green asterisk), and the mass hierarchy is normal. As in the case shown in Fig. [71 allowed 
regions exist also in the inverted mass hierarchy regime when only a single detector at L = 3000 
km is considered, see the panels in the second row. 



similar analysis to understand the feature of the degenerate solutions by using the analytic 
framework and the results are presented in Appendix O 



VI. DISCOVERY POTENTIAL OF NSI PARAMETERS 

In this and the next sections, we discuss the discovery potential of the effects induced 
by NSI, and by the SI, respectively. Throughout these sections we present the sensitivity 
contours which are calculated with two different choices of the assumptions. Choice R and 
Choice O defined in Sec. IIIII We compare the results of the two choices of the experimental 
uncertainties. Typically, the sensitivities obtained with Choice R are worse by a factor of 
several to an order of magnitude than the ones with Choice O. 
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FIG. 12: An example of Type B. Allowed regions with. Cer- 

The input parameters are taken as: 



sin^ 29,-, = 10" 



2 X 10 , (5 = 37r/2, and (per = 77r/4 and normal mass hierarchy. 
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Example of Energy Spectra of Same Sign Degenerate Solutions 
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Example of Energy Spectra of Mixed Sign Degenerate Solutions 
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FIG. 13: Energy spectra of the oscillation probability for a system with non-zero NSI element £e/^ 
corresponding to the same Am|]^-sign intrinsic degeneracy (left panel) and the flipped Am|^-sign 
degeneracy (right panel). 



A. Discovery potential of NSI 

Let us discuss first the sensitivity to discover the non-zero effect of NSI. We determine, 
for some fixed values of ^13 and S, the 2 (3) a CL regions of the discovery potential of NSI 
by the condition, Xmin(^ = 0) ~ Xminl^^^^^ value of e and 0) > 4(9), for 1 DOF, by varying 
freely ^13 and 6 and the choice of mass hierarchy in fitting the input data. If we cannot fit 
the input data generated with non- vanishing Sea by hypotheses with arbitrary values of the 
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SI parameters without NSI, we conclude that we can disentangle NSI effects from those of 
standard oscillations. Note that the region determined by this condition includes the case 
where is totally undetermined or unconstrained. 



Discovery potential of for sin 29j-| = 10 and5 = Ji Discovery potential of for sin' 




FIG. 14: Regions where the non-zero NSI effect caused by £e^ can be identified for the case 
sin^ 2^13 = 10~^ and 6 = it (left panel) and 6 = 37r/2 (right panel). The upper and the lower 
panels are for the cases with measurement at the near (L = 3000 km) detector alone and the one 
combined with the far (L = 7000 km) detector, respectively. The red solid and the blue dotted lines 
are for the sensitivities at 2a and 3a CL, respectively. The thick and the thin lines are for Choice 
R (with systematic uncertainties and efficiencies) and Choice O of the analysis, respectively, 
defined in Sec. Illli 



Discovery potential of E^^ for sin"20j^ = 10 and 5 = 7i Discovery potential of for sin' '29j3 = 10 and 8 = 3jt/2 
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FIG. 15: The same as in Fig. [H but with sin^ 26*13 = 10~^. 

In Figs, [m and [15] we show the regions where the non-zero NSI effect can be discovered in 
the plane of true values of and \ee^\ for the case where sin^ 26*13 = 10~^ and sin^ 26*13 = 
10~^, respectively, for 6 = n (left panels) and 5 = 37r/2 (right panels). The same values of 5 
are consistently used throughout this section for presentation of the allowed contours. These 
two values are chosen because one is CP conserving and the other CP violating. We think 
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that the choice is a rather conservative one, avoiding the point with the best sensitivities. 
In both figures we have taken the input mass hierarchy to be the normal one. 

First, we note that the sensitivity to \ee^\ by the near detector at 3000 km is good enough 
so that adding the contribution of the second detector at 7000 km does not improve much the 
sensitivity to Se^. The feature is indeed expected as discussed in Sees. lITlandHVl Second, the 
sensitivities to \Be^\ depend very much on (p^fi, and varies from ~ 10"'' to ~ 10~^. We have 
checked that the same general behavior is obtained for the case of inverted mass hierarchy 
as input. We note that in the case where detector efficiency, background and systematic 
uncertainties are considered, indicated by the thick curves in Fig. [TH the sensitivities are 
worse by a factor of a few to one order of magnitude compared to the case where they are 
ignored. 



Discovery potential of e^^ for sin 20^^ = 10 and 6 = 71 Discovery potential of e^^ for sin '26^, = 10 and8 = 3it/2 




FIG. 16: Regions where the non-zero NSI effect caused by Ser can be identified for the case 
sin^ 2^13 = 10~^ and 6 = tt (left panel) and S = 3tt/2 (right panel). The upper and the lower 
panels are for the cases with measurement at the near (L = 3000 km) and the far (L = 7000 km) 
detectors, respectively. The meaning of the color-type of the lines are the same as in Fig. [T31 
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In Figs. [16] and Fig. [T7]we show similar plots but for the system with Ser^ ■ In this system 
the sensitivity for the detector at L = 3000 km has a rather complicated structure which 
depends nontrivially on 0er and 5. It is clearly very difficult to discover l^erl ^ few 10~^ with 
the near detector alone. However, as expected, the synergy of combining the two detectors 
is so strong that the final sensitivity is comparable to the case of Sf,^^ i.e. I^erl ^ a few 10^^. 
The sensitivity to also depends on the true value of ^e,-, but the dependence is weaker 
than that for ^g^. Again, we have checked that the same general features are obtained for 
the discovery potential of e^r in the case of inverted mass hierarchy. Here again the role 
of detector efficiencies, background and systematic uncertainties at the level exemplified in 
this paper, is to worse the sensitivities several times compared to the case where they are 
ignored. 

We note that the results shown in this section can be compared with the ones obtained 



in Refs. |60|, l6ll where the discovery potential of NSI effect by a future neutrino factory was 
studied. The authors of Ref. [6l| considered both e^p, and e^t (but only one of them was 
considered at a time) but for the case where only one detector at L = 3000 km is assumed, 
and it is found that for sin^ 26'i3 = 10~^, the sensitivity to ^g^i and Sf-r are, respectively, a few 
X 10~^ and a few x 10~^ under a similar experimental setup to ours. By taking into account 
some differences of the analysis procedure and assumptions, results of our analysis for the 
detector at 3000 km alone seem to be consistent with the ones found in this reference. 

On the other hand, the authors of Ref. [60| studied how one can optimize the neutrino 
factory setting which gives better sensitivities to constrain NSI parameters assuming two 
detectors with different baselines. In this reference, the sensitivity study was done by varying 



baselines to the two detectors and it was found that the standard set up in j55| of two 
detectors at 4000 km and 7500 km, which is similar to the set up considered in our previous 
paper [l9|] (and employed here), works quite well for determining (or constraining) both 
standard and non-standard neutrino parameters. The sensitivities to \Ee^ and leerl found in 
Table 1 of [gO] which are, respectively, ~ 6 x 10^'^ and ~ 2 x 10~^ seem to be consistent with 
ours if we take into account some differences between their and our assumptions (higher 
efficiencies, higher muon energy, larger total number of useful muon decays, etc.) in this 
work. 



B. Discovery potential of the non-standard CP violation 

From the results shown in the previous section, we have some ideas about the pa- 
rameter regions where the non-zero NSI effect can be established. Let us now ask an- 
other question; For which values of the NSI parameters \e\ and 0, can we establish a 
new type of CP violation due to NSI? Hereafter, we denote such CP violation as the 
non-standard CP violation. To establish the discovery potential for such a new CP vi- 
olation effect, we determine the region of sensitivity with 2 (3) o CL by the condition. 



^ The readers may wonder if the jagged behavior in some of the panels in Figs. [16] and [T7| might be due to 
the artifact of a too crude mesh of the parameter space used in the numerical analysis. We have checked 
that the small island structure does not go away, though with modified shape or merging, even if we 
increase the number of grids used in scanning the parameter space by a factor of 10 to 15. The figures 
presented are the outcome of such an improved analysis. Therefore, we believe that some of the jagged 
behavior are real though the possibility of it being a numerical artifact cannot be completely ruled out. 
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Xmin(0 = or vr) — x^jj^(true value of e and 0) > 4(9), for 1 DOF, by varying freely all the 
parameters except for in fitting the input data. 

In Fig. [18] (Fig. [T9l) we show the region where non-standard CP violation associated 
with (bf,,, can be established for sin 2^13 = 10-3 (sin2 2^i3 = 10"^) and 5 = vr (shown 
in the left panel), or 5 = 37r/2 (in the right panel). We observe that in addition to the 
significant improvement of the sensitivity in the case where two detectors are combined (see 
lower panels), there are qualitative differences in the behavior of the sensitivity contours for 
L = 3000 km and the combined case. 

Let us see the right panel of Fig. [181 the case without systematic uncertainties. At 
L = 3000 km, in spite of the maximal non-standard CP violation {(pe/j. = 7r/2, 37r/2) we 
cannot establish the non-standard CP violation even though {Senl is large. Though it might 
look curious it is easy to understand why this feature arises; It is due to the 0-degeneracy 
(see Sec. IV7^ . Using (^^) with input 6 = 37r/2, sin^ 2^13 = lO'^, and E = 20 GeV, one 
can show that the clone solution of 0i = 7r/2 (37r/2) is located very close to 02 = (vr). 
So we are unable to distinguish the maximal CP violation and CP conservation owing to 
the 0-degeneracy. Fortunately, the degeneracy is resolved when informations from the far 
detector is combined. 



Discovery Potential of Non-Standard CPV due to (|)^ Discovery Potential of Non-Standard CPV due to (|)^ 




FIG. 18: Regions where the non-standard CP violation caused by 0e^i 7^ or 7^ vr can be 
established for the case sin 29i3 = 10-3, 6 = TT (left panel) and (5 = 3vr/2 (right panel). Mass 
hierarchy was taken to be normal in the input. The meaning of the color- type of the lines are the 
same as in Fig. [T31 



After combining the results from two detectors, we conclude that if 0.2 < (pefi/T^ ^ 0.8 or 
1.2 < 0e/i/7r ^ 1-9 we can identify the effect of non-standard CP violation down to l^g^l ~ a 
few to several x 10"^ at 3 cr CL for Choice R of the analysis, depending on the values of S 
and sin^ 26*13. For Choice O the sensitivity is better by a factor up to several depending upon 
the SI parameters. We observe that if ^e^i/^r ^ 0.2 or 0.8 < (pe^/T^ ^ 1-2 or ^e/x/^r ^ 1-9 it 
seems practically impossible to establish non-standard CP violation. Here the input choice 
for the mass hierarchy is normal, but the results are qualitatively the same in case we take 
the input mass hierarchy to be inverted and simultaneously change the input 5 — vr — 5. 

In Figs. [20] and [21] we show similar plots but for the discovery of non-standard CP 
violation associated with 0er- We note that for these cases the effect of the synergy of 
combining two detectors is even larger, as expected. We observe that, compared to the 
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FIG. 19: The same as in Fig. [18] but for sin^ 26113 = 10" 



Discovery Potential of Non-Standard CPV due to 



10 ' - 



2 a (w. sys) 

3 a (w. sys) 

2 CT (w/o sys) 

3 CT (w/o sys) 



Normal (input) 



10' 



^ 3 10 




Discovery Potential of Non-Standard CPV due to 



2ct (w. sys) 
3 CT (w. sys) 



Normal (input) 



sin"29, = 10 \ 6 = 3jt/2 - 



FIG. 20: 

0er / TT. 



The same pfots as in Fig. [18] but for the non-standard CP violation caused by (per 7^ or 



system, despite the discovery potential of non-standard CP violation with the detector at 
3000 km alone is rather poor especially for 0er/7r ^ 1, after combining two detectors, the final 
sensitivities are not very much different from that for the system. For 0.1 < (Pst/tt < 0.8 
or 1.2 < (per/'^ ^ 1-9, with Choice R, we can identify the effect of non-standard CP violation 
down to Iserl ~ 10"^ or smaller at 3 cr CL, depending on the values of 6 and sin^ 26'i3. Similar 
to the case for the system with ^e^, for ^et/tt < 0.1, 0.8 < ^et/tt < 1.2, or ^et/tt > 1.9, it 
seems practically impossible to establish non-standard CP violation. 

The results for the sensitivity to non-standard CP violation for the inverted mass hier- 
archy is similar to that for the case of the normal mass hierarchy shown in this paper. 

We observe that including the detector efficiencies, event backgrounds and the systematic 
uncertainties reduce the sensitivities to NSI and the non-standard CP violation by a factor 
of a few to several (even an order of magnitude in a rare case), as expected. But, their 
effects do not appear to change the qualitative features of the sensitivity contours. 

Our results for Eer can be compared with the ones found in j20|. From Fig. 1 of this 
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FIG. 21: The same as in Fig. [20] but for sin^ 26*13 = 10~^. 



reference the non-standard CP violation (due to 0er 7^ 0, vr) can be estabhshed at 3 cr for 
l^erl larger than ~ (7— 10) x 10^^ when the true value of is not so close to or vr for the 
case where the true value of 6*13 
with ours with Choice R. 



0. The sensitivity obtained in [20] is in good agreement 



VII. IMPACT OF NSI ON STANDARD CP VIOLATION AND MASS HIERAR- 
CHY 

At this point it is important to examine whether NSI can obscure the discovery of the 
standard CP violation due to 6 and the neutrino mass hierarchy, and if yes to what extent. 



A. Impact of NSI on the establishment of the standard CP violation 

First let us consider the sensitivity of our setup to 6 without NSI. We show, in Fig. 
in the plane of true values (input) of 6 and sin^ 2^13, the regions where CP violation can be 
established at 2 and 3 a CL for 1 DOF. For this plot we assumed that there is no effect 
of the NSI both in the input data and in the fit. We determine the 2 (3) a CL regions by 
freely varying ^13, 6 and the mass hierarchy in fitting the input data with the condition, 
Xmin(^ = or vr) — x^jj^(true value of 6) > 4(9), for 1 DOF. The upper panel in Fig. [22] 
shows the case where only the detector at 3000 km is considered whereas the lower panel 
is the case corresponding to the combination of detectors at two baselines. Note that, if 
6 is very close to or tt, it is impossible to establish CP violation, a well known fact. In 
the case where only one detector at 3000 km is assumed, there is a small region at around 
6 ~ 37r/2 and sin^ 2^^i3 ~ 3 x 10~^ where the sensitivity is significantly reduced. The loss 
of CP sensitivity at these particular values of 6 and ^13 occurs when the mass hierarchy is 
unknown. This fact has been noted before, see Fig. 85 in Ref. jHI and comments in this 
reference. 

In passing we remark that there exist significant differences between the sensitivity regions 



given in upper panels of Fig. [22] (i.e., blue dotted 3a curves) and Fig. 85 in Ref. [55||. We note. 



however, that if we assume the same detection efficiency as used in Ref. j55|, we obtained 



28 



<x> 

(N 



10^ 



10" 



10 



10"- 



10 



0.0 



Standard CPV sensitivity without NSI 



T 



2 a (w 


sys) 


3 o (w 


sys) 


2a(w 


sys) 


3 a (w 


sys) 


3 a (w 


sys, low cff) 



^ — HH- 



H — f-H 




-T ^ , , r- 
Nornial (input) 



PH — h 



— h 



0.5 



1.0 

8/71 



1.5 



FIG. 22: Sensitivity to discovery of standard CP violation. Here no effect of NSI is assumed both 
in the input data and in the fit. The upper panel shows the case where only the detector at 3000 
km is considered, whereas the lower panel is the case corresponding to the combination of detectors 
at two different baselines. For the sake of comparison, only for the case of 3000 km, we also show, 
by the green dashed curves, the case with lower detection efficiencies as for Fig. 81 of 55l |. 



Standard CPV sensitivity with non-zero £. (in the fit), e, (input)=0 



Standard CPV sensitivity with non-zero £ (in the fit), (input)=0 



10"' 



.9 



2 CJ (w. sys) 

■ ■ ■ ■ 3 a (w. sys) 

2 a (w/o sys) 

3 a (w/o sys) 




FIG. 23: Similar plots as in Fig. [22] but with non-zero NSI allowed in the fit; The input data was 
generated without NSI but non-zero values of (left panel) and Eer (right panel) were allowed 
in the fit. 



roughly the same regions as shown by the green dashed curves in Fig. [221 In comparing 
our results to the one in [551] we should also take into account the difference in baseline; 
Ours is 3000 km whereas the one in Fig. 81 of Ref. 55| is 4000 km. 



As can be seen in Fig. 81 of ^55|], the efficiency used in this reference is roughly a factor of 2 lower than 
the one we take in this paper. Our efficiency is an energy independent approximation to the one given in 



Ref. 



56| 
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Let us now discuss how this result can be affected by the presence of NSI effects below 
the discovery sensitivity. We show in Fig. [231 similar plots as seen in Fig. [22] but assuming 
the presence of non-zero Sen (left panel) and Est (right panel) in the fit. At each point of S 
and sin^ 26*13 in this plot, as in the previous case, we have generated the input data without 
NSI but in fitting the data, we vary freely not only 6, ^13 and the mass hierarchy choice, but 
also the values of Sgfj, (or of egr)- This is done in order to see to what extent the presence of 
NSI in the fit aggravates the discovery potential for (standard) CP violation. 

For the case where non-zero Eefj, is assumed in the fit, by comparing the results shown in 
Fig. [22] with the ones shown in the left panel of Fig. [231 we can conclude that the sensitivity 
regions for discovery of standard CP violation are not strongly affected by the presence of 
small NSI effects (below their discovery level). Neither if we have only one detector at 3000 
km but particularly so if an additional one at 7000 km is combined. However, this conclusion 
is somewhat different when a non-zero Eer is assumed in the fit, as we can see in the right 
panel of Fig. [23 In this case, we conclude that the region of sensitivity to standard CP 
violation is significantly diminished if only the near detector at 3000 km is considered. It 
is even more so for Choice R of the analysis. Fortunately, a similarly good sensitivity 
to CP violation as in the case without NSI is recovered after the combination with the far 
detector, as can be seen in the right-lower panel in Fig. [23] 

We have repeated the same exercise for the case where the input mass hierarchy is in- 
verted. In this case, even in the standard scenario (no NSI at all) the sensitivity region 
is reduced if only the near detector is considered, but if one includes the far detector the 
sensitivity becomes comparable to the one shown in Fig. [25] for the normal hierarchy. In fact 
it is even better for 5 ~ and slightly worse for 5 ~ vr. Here, again the effect of non-zero 
Egfj, and Eer at the single or the combined baselines is similar to the one for the normal mass 
hierarchy. 

Let us also consider the case where the input values of NSI parameters take non-zero 
values which are within the discovery reach. Here, as in the previous case, we vary freely 
the SI and NSI parameters in the fit. We have obtained qualitatively similar results to those 
without input value of NSI (since they are similar, we do not show the plots). We observe 
in this exercise that the difference from the previous results without input NSI at 3000 km 
is quite small for E^fj,, but it is not so small for the case with Eer for the values of the input 
NSI parameters considered. Fortunately, when two detectors are combined, the impact of 
non-zero NSI on CP violation sensitivity does not appear to be so large within values in the 
discovery reach of our setup. We have verified that for the inverted mass hierarchy in the 
case on non-zero input NSI we also reach very similar conclusions. 

B. Impact of NSI on the resolution of neutrino mass hierarchy 

It is generally believed that the neutrino mass hierarchy can be determined "relatively 
easily" in neutrino factory measurements because of sufficiently strong earth matter effect 
due to the long baseline. This is in fact true when only SI are at play. However, once a new 
type of degeneracy, which involves both the normal and the inverted mass hierarchies, is 
uncovered in systems subject to NSI, it is legitimate to ask whether it can affect resolution of 
the mass hierarchy, and if yes, how seriously. Because of the feature discussed in Sec. IIVD[ 
we anticipate that the presence of Egr can be harmful. Even if we assume l^erl to be below 
its discovery potential value, the solution of the mass hierarchy with a single detector at 
3000 km is unsure and the need for a second detector at 7000 km seems to be imperative. 
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We study the region in the input parameter plane 5 — sin^ 2^13 where the neutrino mass hi- 
erarchy can be determined by generating a data set assuming a certain input mass hierarchy 
and trying to fit the data with opposite hierarchy. We say that the mass hierarchy can be de- 
termined at 2 (3) 0" CL for 1 DOF if Xmin(oPPOsite hierarchy) —Xmin(i^Pi^t hierarchy) > 4(9). 

First let us investigate the sensitivity without considering NSI, which means that it is 
absent not only in the input but also in the fit. In Fig. [21] we show the sensitivity region in 
the (5— sin^ 2^13 plane where the mass hierarchy can be determined at 2 and 3 a CL for 1 DOF 
with a single detector at 3000 km. In the left and the right panels the input mass hierarchy 
was taken to be normal and inverted, respectively. We observe that in both cases there are 
small regions in this plane in which the mass hierarchy cannot be determined. However, we 
have checked that if an additional detector at 7000 km is included in the analysis the mass 
hierarchy can be solved in the entire plane covered in the analysis. For the mass hierarchy 
determination, the impact of the non-perfect detection efficiencies, non-zero backgrounds 
and systematic uncertainties is larger compared to the standard CP violation sensitivities 
(see Fig. 
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FIG. 24: Regions in the 5 versus sin^ 2^13 plane where the neutrino mass hierarchy can be estab- 
lished at 2 and 3 a CL for 1 DOF with a single detector at 3000 km, for the case where no NSI effect 
is present. In the left (right) panel the input mass hierarchy is taken to be normal (inverted). It 
can be shown that the mass hierarchy is resolved in the whole region if the far detector is combined. 



Let us now switch on NSI effects and see their impact on the mass hierarchy determina- 
tion. In Fig. [25|, we show similar plots as in Fig. [2l]for non-zero Sen- In the upper panels of 
Fig. [25] we have considered the input to be below its discovery limit (null for all practical 
purposes) but allowed non- vanishing e^fj, in the fit. In the lower panels of Fig. [2S]we take the 
input to be in the discovery region of our setup, l^e^l = 10"'^ and (p^^ = 77r/4, and vary 
freely the SI and NSI parameters in the fit. In both cases, in the region of £e^j below or within 
discovery reach, we observe a slight but not very significant decrease of the sensitivity for 
both the normal (left panels) and the inverted (right panels) mass hierarchies. The regions 
where sensitivity loss is observed are limited to the area around no sensitivity region without 
NSI in Fig. [211 We note that the impact of turning on the systematic uncertainties in the 
manner of Choice R produces a sizable effect. We note that if sin^ 2^13 < 10~^, for some 
range of S the mass hierarchy can not be determined even if we include the far detector at 
7000 km (see the bottom panels of Fig. [25]) . 

Now let us examine the case where is turned on, whose results are shown in Fig. [221 
The organization of Fig. [2^] is as follows; The 1st (2nd) row panels are for case where no 
NSI, or the one below discovery limit, is assumed in the input, for the near detector alone 
(near and far ones combined). The 3rd (4th) row panels are for non- vanishing input value 
of Eer within the discovery region both at the near detector alone (combined ones). 
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0.0 0.5 1.0 1.5 2.00.0 0.5 1.0 1.5 2.0 

5/71 8/71 



FIG. 25: Same as in Fig. [24] but for non-zero £e^i- In the upper panels the input Sefj, = 0, but we 
allow it to be non-zero in the fit. In the lower panels the input is \£e^l\ = 10~^ and 0e/i = 77r/4. 

Let us first look at the case where no NSI is considered in the input. If we consider only 
the near detector at L = 3000 km, the region where the mass hierarchy can be disentangled 
severely shrunk when NSI is switched on. It is the case for both input hierarchies. We 
note that, compared to the case of Se^, the effects of the systematic uncertainties are least 
prominent for the near detector alone (see the 1st row of Fig. [26l) . It becomes impossible to 
resolve it if sin^ 2^13 < 10^'^, and the sensitivity is lost for some values of 5 at larger values 
of sin^2^i3. However, by looking the 2nd row of Fig. [22] we observe that the power of the 
synergy of the two detectors is quite strong so that the combination of the two detector 
allows us to determine the mass hierarchy for most of the parameter regions we considered 
in this paper even in the presence of very small Eer- 

Then, what happens if (input vale of) Eer is not so small? As seen in the 3rd row panels 
of Fig. [26], it strongly aggravates the mass hierarchy determination. For instance, if the 
input hierarchy is normal (left panel) the hierarchy is undetermined if sin^ 26'i3 < a few 
xlO^^. If the input hierarchy is inverted (right panel) sensitivity to the hierarchy very much 
depends on 5 and it can be determined down to sin^ 2^13 < 10^^ only in a limited region of 
5. In the bottom panels of Fig. [22] we show the result of combining the near and the far 
detectors. Even with this combination there is a small region where the mass hierarchy is 
not determined. 

We note that compared to the Ee^ element, the impact of the detection efficiencies, back- 
ground and systematic uncertainties is relatively small. It appears that this is partly because 
the impact of non-zero Eer itself is already large before inclusion of the uncertainties. 
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Sensitivity to Mass Hierarchy with e NSI effect 
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FIG. 26: Similar plots as in Fig. [25] but for non-zero SeT- In the upper panels the input e^t = 0, 
but we allow it to be non-zero in the fit. In the middle panels the input is legrl = 2 x 10~^ and 
(/>eT = 7r/4. In the lower panels we present the result of the combination of the detectors at 3000 
and 7000 km for the input (|eer| =2 x 10"^, c^er = vr/4). 



VIII. CONCLUSIONS 

In this paper, we have studied the question of how to distinguish between physics effects 
due to SI and NSI in neutrino oscillations. They include, most notably, discriminating 
effects of possible non-standard CP violation due to phases associated with the NSI elements 
from the standard CP violation caused by the lepton KM phase 5, and vice versa. They 
also include a related question of how NSI affects the determination of the neutrino mass 
hierarchy. Our study was done in the context of a future neutrino factory endowed with an 
intense muon storage ring of 50 GeV, delivering 10^^ useful /x-decays a year and operating 
for 4 years in neutrino and 4 years in antineutrino modes. We assumed two magnetized iron 
detectors, one at 3000 km and the other at 7000 km from the neutrino source, each with 50 
kton fiducial mass. We considered only data from the golden channels v^. —>■ and Vf. ^ t^h- 

In this work we have considered NSI effects only in the neutrino propagation, and turned 
on only one of the relevant NSI elements for the golden channel, Se^ = \ee^i\e^'^''^ or Ser = 
l^erle*'^''^. Since neutrino oscillation in such a system is very complicated we have studied 
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the effect of a single e (two real parameters) at a time. We have fixed the values of the 
standard oscillation parameters at their current best fit values, expect for sin^26'i3, 5 and 
sign of Am|^, which we also vary freely in our analysis. 

Prior to the presentation of the full discovery potential of NSI and SI parameters we have 
analyzed the structure of neutrino oscillation with NSI. We have utilized the bi-probability 
plot in P(z/e — > Vfj) — P{i'e ~^ i^fj.) space, which are drawn by varying either NSI phase 
or the KM phase 6, to illuminate the characteristic properties of the neutrino oscillation in 
systems with NSI. It revealed to be a powerful tool to understand the synergy between the 
two detectors, and the difference between the Eefj. and Eer systems. Furthermore, we have 
given an overview of the sensitivities to NSI and the SI parameters by classifying the data 
set of allowed region contours into four different types depending upon the degree of synergy 
between the near and far detector measurements. 

Most notably, we have observed in the near detector measurement the phenomenon of 
parameter degeneracy which is similar to the one in standard neutrino oscillation but with 
active participation of NSI effects. Though it is a highly complicated problem, we were able 
to control it analytically in the restricted setting in the presence of only a single type of 
e. The NSI enriched parameter degeneracy discussed in this paper contains the Am|]^-sign 
flipped degeneracy, the intrinsic one, and the one called the 0— degeneracy. The last one is 
the special case of the more generic cases, but it is best characterized as the 0— degeneracy 
because the solutions are different essentially only by the values of 0. We have shown that 
these degeneracies are very robust against spectrum analysis and so can potentially disrupt 
the resolution of CP violation and the mass hierarchy. However, it is also observed that in 
many cases the parameter degeneracies can be resolved by adding the far detector. This is 
why, especially for NSI due to Eer, the inclusion of the 7000 km detector is imperative. 

We have analyzed the discovery potential of NSI with our setup, by performing two typical 
cases of analyses. Choice R and Choice O as defined in Sec.lIIIl with and without detection 
efficiency function, background and systematic uncertainties, respectively. In the range of SI 
parameters covered in this paper and independent of the neutrino mass hierarchy, we have 
established that measurements performed by a single detector at 3000 km can discover NSI 
due to Ee^ down to l^e/il ~ 10~^ — 10""^ (Choice O) and a few to several times larger than 
these for Choice R, the exact value depending on 0e^. In the e^^ system, inclusion of the far 
detector does not improve this much. However, in the case of the discovery of NSI due to E^r 
the synergy between the two detectors is very strong. Although one can only hope to reach 
down to {EetI ~ a few xlO~^ (Choice R) with the near detector at 3000 km, the inclusion 
of the far detector makes the sensitivity to Eet similar to that for Eefi- We have identified 
through the discussions in Sees. [TTl and [IV] why such disparity between sensitivities in the 
Eefj, and the Eer systems arises. We note that taking into account the non-perfect detection 
efficiencies, non-zero backgrounds and systematic uncertainties change sensitivities a few to 
one order of magnitude. 

We have also investigated the potential to unravel CP violation associated with NSI 
using such an experiment. We have concluded that to establish CP violation associated 
with 0e^ the addition of the second detector indeed helps. For both mass hierarchies, if 
0.2 < 0e^/vr < 0.8 or 1.2 < 0e/x/7r < 1.9 we can discover non-standard CP violation pro- 
moted by this phase down to \Eefj.\ ~ (a few - several) x 10^^ (for Choice R, depending on 
sin^ 29i3 and 5) at 3 cr CL. For CP violation associated with 0er the synergy between the two 
detectors is more efficient than the one in the (pefi system. The near detector alone gives poor 
sensitivities but after combining two detectors, if 0.1 < 0er/7r ^ 0.8 or 1.2 < 0er/7r < 1.9, 
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we can identify the effect of non-standard CP violation down to |£er| ~ a few x 10~^ — 10~^ 
(for Choice R, depending on sin^ 2^13 and 5) at 3 a CL, independently of the neutrino 
mass hierarchy. Improvement of the systematic uncertainties does affect the sensitivity to 
non-standard CP violation by a factor of a few but not so dramatically. 

We have also checked to what extent the existence of NSI can aggravate the discovery 
potential of standard CP violation in the S — sin^ 2^13 plane. We have observed that NSI 
effects induced by Sg/* (either with a magnitude below or above the discovery reach of our 
setup) will not spoil much the discovery potential of standard CP violation with respect to 
the standard case, even if we have a single detector at 3000 km. Nevertheless, NSI effects 
induced by Ser are potentially harmful if only the near detector is considered, even for values 
of Iserl below the sensitivity reach. Fortunately, the discovery potential possessed by the 
setting without NSI is essentially recovered by the inclusion of an additional detector at 
7000 km. We have checked that these conclusions do not depend on the input neutrino mass 
hierarchy. 

Finally, we have studied the impact of NSI on the resolution of the neutrino mass hierar- 
chy. For non-zero 8^/^, with the near detector at 3000 km, we observe a slight but not very 
significant decrease of the region in the 5 — sin^ 2^13 plane where the neutrino mass hierarchy 
can be established. By adding the far detector at 7000 km, one can distinguish the mass 
hierarchy for almost all values of S if sin^2^i3 > 10~^. For non-zero Eer, the region where 
the mass hierarchy can be determined with the near detector alone severely shrinks. Here 
again, the power of the combination of two detectors is sufficiently strong, allowing the mass 
hierarchy to be determined in almost the whole parameter space of 5 and ^13 considered in 
this work (except for a small region if Iserl is rather large). We observe that the impact of 
non-perfect detection, background and systematic uncertainties on the determination of the 
mass hierarchy is somewhat larger for the case where ^e/* is present than for the case where 
Eer is, though the effect of the latter is already large without these factors. 
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APPENDIX A: EXPRESSION OF THE APPEARANCE OSCILLATION PROB- 
ABILITY FOR THE BI-PROBABILITY PLOT 



The appearance oscillation probability in the neutrino channel in system with 
/i, r) are given by 
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The analogous expression of P(z/e ^ u^) with eg,- can be obtained from ([3]), or by noting 
that they originally appear in the particular combination of the generalized atmospheric and 
the solar variables ^4| though it is obscured a little in ([3]). (See Eq. (6.5) in [33].) The prob- 
ability for anti-neutrino can be obtained from the one for neutrinos by doing transformations 
a —a, 5 —5, and 6 —d). 
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APPENDIX B: UNDERSTANDING ROTATING ELLIPSES IN THE BI- 
PROBABILITY PLOT 



It can be easily shown that the bi-probabihty trajectory drawn by varying holding S 
fixed (or, vice versa) takes the form of ellipse. We obtain the equation which determines the 
major (minor) axis from (jl]) 
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where 



S = TZ + B cos S + CsinS 
T = J — Bsm6 + C cos 6 
S = TZ + B cos 6 — Csm6 
T = i + Bsm6 + C cos6. 

Hence, the slope of the major axis, Ovaryc/., is given by 
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so the slope of the axis is independent of the value of l^eal- 

The 6 dependence of ctvaryt/) are plotted in Fig. [271 The slope of the major axis is mostly 
negative and positive in the and Eer systems, respectively. When 9i3 is small, the 
coefficients of 6 independent term [TZ and X) relatively play the dominant role in fIBSI) . the 
rotating behavior of the ellipse is soft. On the other hand, the slope of the major axis 
changes dramatically in the region around 6 = it with large 6'i3 as we saw in the Figs. [1] and 
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FIG. 27: 



S S 
The slope of the major axis of the ellipse which is made by varying 



in the ee/i(ee 



system, sin 2^13 = 10 (10 ) for the red soUd (blue dashed) line, ayary^ is independent of |ee 



Similarly, we can work out the behavior of slopes as a function of in the case of the 
bi-probability trajectory drawn by varying 6 by holding 0. Skipping details, we present the 
results of the slope ftvaryS of 5-varied bi-probability plots in Fig. [281 We observe the similar 
rotating behavior of the bi-probability diagram as in the (5-varied ellipses. 
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FIG. 28: The slope of the major axis of the elhpse which is made by varying 6 in the ee/x(£er) 
system, leg^tl = 5 x 10^'^(5 x 10""^) and |eer| = 2 x 10~^(5 x 10~^) for the red sohd (blue dashed) 
line. Ovary (5 is independent of 613. 

APPENDIX C: VARIOUS DEGENERATE SOLUTIONS AND ANALYTIC SO- 
LUTIONS 

We note that for a given true solution, it is possible to obtain degenerate solutions 
using analytic expressions as follows. To do this let us take the analytic expression of the 
appearance oscillation probabilities P(z/e —>■ v^) in ([3]) (and its anti-neutrino counterpart) 
with the notation P(z/e — > z/^) = PeAi('5, 6*13, l^eal, 0ea; -E) (a = /i, r). Since the number of 
unknown parameters are four, 6*13, 5, and 0, we need four observable quantities. We take 
the oscillation probabilities P(z/e — >■ v^) and P{i'e i^n) at two different energies Ei and E2 
for these four inputs. Taking the assumed input values of the four parameters given in the 
caption of Fig. [H] and Fig. [121 we solve the equations 

P,,(r-, err, iCl^rr, E,) = ^e,(5^, ^^S, E,) {a = ^^,r■^ = 1, 2) (Cl) 

numerically to obtain the degenerate solutions to where the superscript "D" is attached. To 
solve (ICip we arbitrarily take two reference energies as Ei = 10 GeV and E2 = 20 GeV. 
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solution of (ICTl) 
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approximate solution of (jCip 


inverted 


0.0035 
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0.0046 
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approximate solution of (jCip 


inverted 


0.0025 


1.4 


0.0048 
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TABLE III: Presented are solutions of the degeneracy equation (jCip for input parameters similar 
to the ones used in Fig. [11] but with c/ig^ = 37r/4 as given in the first row. See the text for 
explanation of what "approximate solution of (|Cip " means in the first column of the Table. 

In Tables [Inland [TVlwe present examples of such degenerate solutions for the system with 
Eefj. and Est, respectively. The first column of Tables [TTTI and [IVl is to specify the nature of 
the solutions. The label "approximate solution of fICip " implies the following situation: By 
solving (ICip with the input parameters in the second column we obtain a complex solution 
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nature of the solution 


hierarchy 


sin^ 26*13 


6 


l^erl 


<t>eT 


input (Fig. [12]) 


normal 


0.001 


fvr 


0.02 


h 


solution of (ICip 


inverted 


0.0016 


2.77 


0.022 


0.51 



TABLE IV: Presented are solutions of the degeneracy equation (jCl|) for input parameters corre- 
sponding to Fig. [12] given in the first row. 



nature of the solution 


hierarchy 


sin^ 26*13 


5 






input 


normal 


0.001 





0.005 


1^ = 3.9 


solution of (|C1|1 


normal 


0.001 





0.005 


2.4 


approximate solution of (jCip 


inverted 


0.0034 


1.5 


0.0047 


3.4 


approximate solution of (ICip 


inverted 


0.0025 


1.4 


0.0049 


5.1 


input 


normal 


0.0001 


i7r = 1.6 


0.005 


fvr = 3.9 


solution of ([Cip 


normal 


0.0002 


0.45 


0.004 


1.8 


approximate solution of (1C1|) 


inverted 


0.0015 


0.6 


0.0052 


6.0 



TABLE V: Similar table as TABLE [m] 



which cannot be regarded as the physical one. The solutions given in Table IIIII are real 
numbers which are close enough to the complex solutions. It should also be noticed that 
the degeneracy equations (ICip sometimes have solutions which do not survive in a form of 
allowed contours as a results of analysis of neutrino factory measurement at L = 3000 km. 
For example, there is a solution sin^ 26*13 = 0.0028, 5 = 3.5, \ee^\ = 0.008, and (pefj, = 5.1 for 
the same input as given in Table lllli We have confirmed that this solution indeed solves 
the equation ( ICll) at ii^ = 10 GeV and ii^ = 20 GeV, but the oscillation probabilities deviate 
from the input ones at elsewhere in the energy spectra. Therefore, the degenerate solution 
was lifted by the spectral informations used by our numerical analysis. 

In Tables FVllVIl we also show some other examples. Table |V] is similar table as Table lllli 
Input values of parameters are close to the first solution of Table IIIII Therefor the solutions 
of this set are very similar to Table llllf s. These relationship exists in the system with eer 
too. In Table IVI| input values of parameters of upper set are close to the first solution of 
lower set with similar another solution. 
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nafiirp of fhp ^nlntinn 


Vi iprairnv 


sin^ 26*1 q 


§ 


peT 1 


6 


innn t, 




0.001 


vr = 3.1 


0.02 


iyr = 0.79 


solution of dSIl) 


normal 


0.001 


3.1 


0.022 


5.4 


solution of (ICTI) 


inverted 


0.001 


4.5 


0.016 


3.7 


input 


normal 


0.001 


vr = 3.1 


0.02 


1^ = 5.5 


solution of jCl]) 


normal 


0.001 


3.1 


0.018 


0.63 


solution of dSI]) 


inverted 


0.001 


4.5 


0.016 


3.9 



TABLE VI: Similar table as TABLE [El 
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